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Abstract 

A search for direct production of Higgs bosons in the di-tau decay mode is performed with 
86.3 ± 3.5 pb -1 of data collected with the Collider Detector at Fermilab during the 1994—1995 
data taking period of the Tevatron. We search for events where one tau decays to an electron plus 
neutrinos and the other tau decays hadronically. We perform a counting experiment and set limits 
on the cross section for supersymmetric Higgs boson production where tan (3 is large and m\ is 
small. For a benchmark parameter space point where m^o = 100 GeV/c 2 and tan (3 = 50, we limit 
the production cross section multiplied by the branching ratio to be less than 77.9 pb at the 95% 
confidence level compared to the theoretically predicted value of 11.0 pb. This is the first search 
for Higgs bosons decaying to tau pairs at a hadron collider. 

PACS numbers: 14.80.Cp,13.85.Rm,11.30.Pb,12.60.Fr,12.60.Jv,14.60.Fg 
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I. INTRODUCTION 

The Higgs mechanism in the Minimal Supersymmetric Standard Model (MSSM) 
provides a way to assign a mass to each particle while preserving the gauge invariance of the 
theory just as in the Standard Model (SM). The CP-conserving MSSM contains two SU(2) 
Higgs doublets yielding five physical particles - four CP-even scalars (h°, H°, H~ and H + ) 



lysical 

IBB. 



and one CP-odd scalar (^4°) |4| |5j, |6j] . Here, h° is the lighter of the two neutral scalars. In the 
MSSM there are two parameters at tree level that are conventionally selected to be tan (3 and 
m A o. The parameter tan/3 is the ratio of the vacuum expectation values of the two Higgs 
doublets, and m^o is the mass of the pseudoscalar Higgs particle. When tan (3 is large, the 
cross section for direct production of Higgs bosons through gluon fusion becomes enhanced, 
making that an appealing region for searches at the Tevatron. The coupling strength of 
the A boson to down-type fermions with mass is proportional to mj x tan/3, hence 
couplings to tau leptons are enhanced. The couplings of the ho are similarly enhanced for 
many possible models. 

No fundamental scalar particle has yet been observed in any experiment. The four experi- 
ments at LEP have each performed a search for h°/A° produced in the process: e + e~ — > h°Z 
and e + e~ — > h°A°. The combined results of four experiments have constrained the theory, 
excluding m A o < 91.9 GeV/c 2 , m h o < 91.0 GeV/c 2 and 0.5 < tan/5 < 2.4 at 95% confidence 
level Zt- Another search for h°/A° produced in association with two bottom (b) quarks 
and decaying to two b quarks was performed at CDF earlier 0. This previous search was 
sensitive to the high-tan (3 region, excluding tan/3 > 50 for m A o = 100 GeV/c 2 . 

This paper presents the results of a search for supersymmetric (SUSY) Higgs bosons 
directly produced in proton- antiproton collisions at a center-of-mass energy of 1.8 TeV using 
the 86.3 ±3.5 pb" 1 of data recorded by the Collider Detector at Fermilab (CDF) during the 
1994—1995 data taking period of the Tevatron (Run lb). Although the branching ratio to 
b quarks would be largest, that decay mode would be dominated by QCD background, so 
we search for Higgs bosons that have decayed to two tau (r) leptons. Events are selected 
inclusively by requiring an electron from r — > ev e v T and a hadronically decaying tau (77J 
lepton. This semi-leptonic mode was chosen for this search as a trade-off between the 
distinctive electron signature in a QCD environment and the high branching ratio of hadronic 
tau decays. 
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This is the first time that a search for Higgs bosons has been carried out in the di-tau 



decay mode using data 
taus in the final state 



rom a hadron collider. In Run 1, CDF published other analysis with 



I], llo| . We also demonstrate for the first time from such data the 
feasibility of a technique to reconstruct the full mass of a candidate di-tau system, which is 
only possible when the tau candidates are not back-to-back in the plane transverse to the 
beam. 

The sample that passes the final selection cuts is dominated by Z — > rr events. There 
is no evidence for a signal, so we report a limit on a set of MSSM models in a region of 
parameter space where tan f3 is large because this is where the best sensitivity is achieved. 
Since at the Tevatron most directly produced Higgs bosons would be back-to-back, the 
acceptance is small in the region where the mass reconstruction is a discriminating variable. 
Therefore, limits are reported based on a counting experiment using events from the full 
sample. Then, from a subset of the events where the tau candidates are not back-to-back 
in the transverse plane we extract a mass distribution and perform a binned likelihood to 
demonstrate the capability of that technique. 

II. THE CDF DETECTOR 

This section briefly describes the Run 1 CDF detector with an emphasis on the sub- 
detectors important to this analysis. The CDF detector is described in detail elsewhere jlll 
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CDF used a cylindrical coordinate system with the z axis along the proton beam direction. 
The polar angle (9) was reported with respect to the z axis. Pseudorapidity (77) was defined 
as — ln[tan(0/2)]. Detector pseudorapidity (r] d ) was the same quantity with dependence on 
vertex position removed. The azimuthal angle (<fr) was measured relative to the positive x 
direction. 

The CDF electromagnetic and hadronic calorimeters were arranged in a projective tower 
geometry, as well as charged particle tracking chambers. The tracking chambers were im- 
mersed in a 1.4 T magnetic field oriented along the proton beam direction provided by a 
3 m diameter, 5 m long superconducting solenoid magnet coil. 

In the central region covering \t\\ < 1.1, the electromagnetic (CEM) and hadron (CHA, 
WHA) calorimeters were made of absorber sheets interspersed with scintillator. Plastic 
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light guides brought the light up to two phototubes per EM tower. The towers were con- 
structed in 48 wedges, each consisting of 10 towers in rj by one tower in 0. The mea- 
sured energy resolution for the CEM and CHA were a(E)/E = 13.7%/\/Ey © 2% and 
a{E)/E = 50%/y/E^® 3%, respectively 

The central EM strip chambers (CES) were proportional strip and wire chambers located 
six radiation lengths deep in the CEM (radial distance r = 184.15 cm), where the lateral size 
of the electromagnetic shower was expected to be maximal QQ,0,Q. 



It measured 



the position of electromagnetic showers in the plane perpendicular to the radial direction 
with a resolution of 2 mm in each dimension Q] . In each half of the detector (east and 
west), and for each 15° section in 0, the CES was subdivided into two regions in z, with 128 
cathode strips separated by ~ 2 cm measuring the shower positions along the z direction 
with a gap within 6.2 cm of the z = plane. In each such region, 64 anode wires (ganged 
in pairs) with a 1.45 cm pitch provided a measurement in 0. 

A three-component tracking system measured charged particle trajectories, consisting 
of the silicon vertex detector (SVX') f _the time projection chamber (VTX) and the central 
tracking chamber (CTC). The SVX' |2Qj consisted of four concentric silicon layers sitting at 
radii between 2.36 cm and 7.87 cm and providing r — tracking information only. The VTX 
was positioned just beyond the SVX' in radius and measured the position of the collision 
point along the beam for each event to a resolution of 2 mm. 

Beyond the VTX (radially) was the CTC, a cylindrical drift chamber 3.2 m long in the 
z direction, with its inner (outer) radius at 0.3 (1.3) m. The sense wires were arranged 
into 84 layers divided into 9 "super- layers." Five of the super-layers (axial) contained cells 
with 12 sense wires that ran parallel to the beam and provided measurements in r — 0. 
The remaining four super-layers (stereo) sat between the axial layers in radius, contained 
6 sense wires per cell, and were rotated in the r — z projection by 2.5° with respect to the 
beam to provide measurements in r — z. The transverse momentum resolution of the CTC 
was Spt/pt ^ 0.002 GeV/c _1 x px and when combined with the SVX' tracking information 
when available, the resolution was Spr/pr < 0.001 GeV/cT 1 x p T . 
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FIG. 1: The figure shows the expected mass distribution of Higgs bosons produced directly at the 
Tevatron for m A o = 100 GeV/c 2 and tan/3 = 50. We compare the output of Pythia 6.203 [21 1 
with that obtained from the method described in this paper: folding in a simple relativistic Breit- 
Wigner distribution with the mass-dependent cross section calculation. 

III. MONTE CARLO SIMULATION 



The Pythia 6.203 



2l| event generator is used to simulate signal events and backgrounds 



other than fakes. The Monte Carlo (MC) samples that are generated with the standard 
Pythia package required modifications that are discussed in this section. 

A re-summed calculation of the cross sections for direct Higgs boson production at the 



22 1 , which allows a user to estimate 



Tevatron has been performed using the program HIGLU 
cross sections as a function of mass, tan j3 and other parameters. For a given Higgs boson 
mass HIGLU gives the on-shell cross section only. A SUSY Higgs boson produced at large 
tan f3 has a significant width and a tail at low values of the center of mass energy of the 
parton collision that produced the Higgs boson \[Q 2 (see Fig. [TJ). Therefore, scaling the 
Pythia differential distribution to the HIGLU one would underestimate the cross section 
because the off-shell events would not be accounted for. 

To estimate the total cross section, we first retrieve the on-shell cross section as a function 
of mass in bins 1 GeV/c 2 wide using the HIGLU program. Then, the mass- dependent cross 
section is folded into a relativistic Breit-Wigner shape, with the width proportional to Q 2 , 
from Q 2 = (40 GeV) 2 to Q 2 = (200 GeV) 2 . For m A o = 100 GeV/c 2 at tan/5 = 50, the 
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FIG. 2: The px of directly produced Higgs bosons with mass tua = 100 GeV after a sin A(j) > 0.3 
cut has been imposed. 

MSSM cross-section for A + h° production is 122 pb. 

The rate of Higgs boson production in the region of low \/~Q 2 ~ is a source of significant 
uncertainty for the analysis, particularly at high mass and high tan/5. The low tail seen 
in Fig. ^ originates from a steeply falling cross section folded in with an increase in parton 
luminosities at small momentum transfer, folded in with a broad Higgs boson width. When 
we use the method outlined above to obtain a Q 2 dependent cross section, the size of the tail 
is bounded by those obtained when one generates events at the same parameter space point 
with Pythia and Is ajet Q . We compare the result from the default P ythia output with 
that where we use the method above to estimate the systematic error from these low mass 
tails. At m^o = 100 GeV/c 2 and tan/3 = 50, this uncertainty is 2%. At higher mass and 
higher tan/3, this systematic can become significant. At m^o = 140 GeV/c 2 , tan/3 = 80 the 
systematic is 30%. 

It is important to properly model the boson p? for Higgs boson and Z boson events 
because it impacts the relative rates of back-to-back and non-back-to-back events, and the 
di-tau mass resolution for the latter events. Figure El shows the pt distribution of the A 
(m^o = 100 GeV/c 2 ) after the sin A0 > 0.3 cut is imposed, where A0 is the angle between 
the taus in the transverse plane. We see that this cut is approximately equivalent to a cut 
of pt > 15 GeV/c imposed on the parent Higgs bosons. 

In the high-tan j3 region of parameter space probed in this analysis, the direct-production 
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FIG. 3: Feynman diagram of the direct production process. In the high-tan f3 region of parameter 
space, the gluon fusion process occurs predominantly through a triangular bottom quark loop. 



process occurs predominantly through gluon fusion via a triangular bottom quark loop as 
in Fig. El (the direct production of SM Higgs bosons proceeds predominantly through a t 
quark loop). In the default Pythia, the effect of the lighter quark mass in the fermion loop 
on the pt of Higgs bosons is not taken into account. To correct for this, we use 24J which 
performs a perturbative calculation for the differential cross section da/dpr to order with 
a variable quark mass in the triangular loop. The perturbative calculation is valid in the 
region p^, l " s > 15 GeV/c. As has been pointed out, this cut is nearly equivalent to the 
requirement sinA0 > 0.3. We force agreement between the Higgs boson p? distributions 
in the region p^" s > 15 GeV/c and the result of the program with a 5 GeV/c 2 b quark 
in the fermion loop so that the resulting MC sample will have the proper efficiency for the 
sin A0 cut. To do this, we use a reweighting method known as the acceptance rejection MC 
method 

The pt distributions for the Z bosons from Z/'y* — > e e events has been measured by 
CDF in the region 66 GeV/c 2 < M ee <116 GeV/c 2 Q. Z bosons from Pythia tend 
to have a lower average pt than the measured value. Therefore, we reweight both the 
Z/7* — > t t and Z/'y* — > ee samples to force agreement between the MC events and the 
measured spectrum. Only events that lie in the measured region are subject to rejection. 
This correction makes a significant difference in the relative rates of back-to-back and non- 
back-to-back events. Before the correction, 20% of the Z events from Pythia are in the 
non-back-to-back region (here, defined by > 15 GeV/c). After the correction, 26% of the 
events are in the non-back-to-back region. 

CDF measured the cross section for ZH* — > ee in the same mass window as the t>t distri- 

n 

butions above, and reports this cross section to be 248 ± 11 pb [26J. Assuming universality, 
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we scale the generated Z — > tt sample to this cross section in the measured range. 

Polarization effects, which impact the tau lepton kinematics, are properly taken into 
account with Tauola |2jJ. Simulated Higgs boson events tend to produce one tau with hard 
(high-pr) visible decay products and one with soft visible decay products. Taus produced 
from Z bosons are either both hard or both soft. 



IV. EVENT SELECTION AND EFFICIENCIES 



The search is performed in a data sample of events collected using a high-i^y electron 
trigger. In the following, we first describe the trigger system that an event must pass to 
enter this data sample, followed by a description of the electron and tau identification, and 
finally the event selection and mass reconstruction made at the analysis level. 



A. Trigger 

28j and the energy- 



A three-level trigger system was used to select events for data storage 
dependent efficiency of these triggers has been measured [2^. At Level 1 (LI), the trigger 
requires at least one trigger tower with Et > 8 GeV in the CEM and the efficiency was 
measured to be 100% at the 10~ 3 level. At Level 2 (L2), the trigger requires one calorimeter 
cluster with Et > 16 GeV and the ratio of hadronic energy to electromagnetic energy 
(Ehad/Eem) to be less than 0.125. This cluster must be close in azimuthal angle to a track 
reconstructed in the trigger with px > 12 GeV/c. The triggered sample contains 128,761 
events. The efficiency of the L2 trigger for a good quality electron with Et > 20 GeV is 
91 ± 2%, measured from a Z — > ee control sample. Signal events would pass this trigger 
approximately 20% of the time. 



B. Electron Identification 



able H} are made on the 
we briefly summarize 



To improve the purity of the data, further cuts (listed in 
candidate electron. These cuts are described in more detail in 
them here. We use Z — * ee control samples to quantify the degree of agreement between 
the simulation and the data. 
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E T > 20 GeV 
p T > 13 GeV/c 
E/p < 1.5 
Lshr < 0.2 
Ehad/Eem < 0.05 
| Ax | < 1.5 cm 
|Az| < 3.0 cm 

Xltrip < 10.0 
kvTX — z e\ < 5.0 cm 

kvTxl < 60.0 cm 
Conversion Rejection 
Fiducial cuts on the electron 

TABLE I: Cuts made offline to select the data sample used for the search. The electron identifica- 
tion variables are defined in Sec. \IVBl 

A candidate electron is first identified as a calorimeter cluster in the CEM. An electron 
cluster in the calorimeter is formed by merging seed towers (required to have E? > 3 GeV) 
with neighboring towers in rj (E? > 0.1 GeV required). To be called an EM cluster, E? > 5 
GeV is required. Corrections are made to the energy of an EM cluster to compensate 
for variable response across each tower, tower-to-tower gain variations and time- dependent 
effects. A global correction to the energy scale is also imposed jsij. These corrections are 
typically at the level of a few percent. 

We require an electron candidate to have a CTC track pointing to an EM cluster. The 
highest px track pointing to the cluster is the "electron track," and is required to satisfy 
Pt > 13 GeV/c. The measured direction of the track momentum sets the direction of the 
electron candidate. 

We require that the ratio of the energy deposited in the EM calorimeter to the momentum 
of the electron track is not too large (E/p < 1.5). Also, the lateral profile of the EM shower 
left by the electron tau candidate is required to be consistent with electron shower profiles 
as measured in test beam data (Lshr < 0.2). The ratio of energy measured in the CHA to 
energy measured in the CEM (Ehad/Eem) is expected to be small for an electron from a tau 
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decay; we require Ehad/Eem < 0.05. In addition, the electron track projected to the plane 
of the CES must be close to a CES shower position: \Ax\ < 1.5 cm and \Az\ < 3.0 cm. This 
reduces the background from charged pions produced in the neighborhood of neutral pions. 
We also confirm that the profile of the pulse heights produced by the electron shower across 
CES strips is consistent with electron test beam data: xltrip < 10.0. The electron track is 
required to be consistent with a vertex that lies within 60 cm of the nominal collision point. 
Standard CDF fiducial cuts are made on the electron to ensure that the particle arrived at 
an instrumented region of the calorimeter with good response. We reject electron candidates 
which are consistent with having originated from a photon that converted to an electron- 
positron pair in the detector by removing candidates that leave a low-occupancy track in 
the VTX or that have a nearby opposite-sign track. The opposite-sign track must be within 
90° in <fr, separated in r — <j) from the electron track by no more than 0.3 cm measured at the 
point where the tracks are parallel, satisfying | A cot 9\ < 0.06 where A cot 9 is the difference 
between the values of cot 9 for the two tracks. The electron tau candidate must be isolated 
in the calorimeter and in the tracking system. In the calorimeter, we use the standard CDF 
isolation variable defined by: 

77>cone rjicluster 
-n-iso ^cluster ^ 

where E cone is the energy deposited in the calorimeter in a cone of 
AR = a/ (A0) 2 + (A77) 2 = 0.4 around the electron tau candidate and E clustcI is 
the energy of the EM cluster. We require i?i SO < 0.1. We define Ni SO to be the number of 
tracks with > 1 GeV/c within AR < 0.524 of the EM cluster. This cone size is 30°, 
chosen to be the same as the outer radius of the isolation annulus used for identifying 
hadronic taus, described below. A track must originate within 5 cm of the electron track 
along the beam line to be counted in the isolation cone. We require that Ni so = 0. 

We require at least one EM cluster in the event with E? > 20 GeV passing the electron 
identification requirements just described. We refer to this as the electron tau (r e ) candidate. 
If there is more than one r e candidate in the event, we select the candidate with highest Et- 

We correct for inadequacies in simulation of electron identification variables by applying 
a scale factor of 0.869 ±0.016. The scale factor was determined using Z/'-f* — > ee events and 
we have determined that is not pt dependent. 

Since the vertex position affects the acceptance, we reweight the events to force agreement 
between the distribution of primary vertex positions from the simulation and that measured 
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E T > 10 GeV (jet) 
p T > 10 GeV/c (track) 

N trks = i 3 
- 1 ' cone 

jytrks _ q 
aim " 

Fiducial requirements 
1-^VTX — z Th \ < 5.0 cm 



\z Th \ < 60.0 cm 



iV ccs < 3 

A " cone u 



iV ccs = 

ann w 

m T < 2.0 GeV/c 2 
£ > 0.15 



< 0.1 and 1^ < 0.1 



TABLE II: Cuts made offline for tau identification. These variables are defined in Sec. U V CI 



from the Run lb data sample. 



C. Hadronic Tau Identification 



For taus coming from Higgs boson decays, the hadronic tau decay products will be col- 
limated, with an angular deviation from the direction of the tau parent of no more than 
A</> < m T /E r which is ~ 10° for a typical E T ~ 10 GeV. In nearly all cases, tau decay 
products will include 1 or 3 charged tracks and < 2 neutral pions, each decaying to two 
photons (all other decay modes have branching fractions of less than half of a percent). 

The cuts used to select a hadronic tau are listed in Table HU The Th identification cuts 
used here are based on those outlined in (3^ . The main differences are noted in what follows. 
The search for a Th begins with identifying a stiff track associated with a jet cluster. We 
require a track with > 10 GeV/c within AR < 0.4 of a jet cluster with > 10 GeV. 
The calorimeter cluster size is AR = 0.4. The track with the highest pt satisfying this 
requirement is called the tau seed. The Et cut is approximately 75% efficient for signal. 
The pt cut is approximately 65% efficient for signal while rejecting 80% of QCD jets (after 
the Et cut has already been imposed). 
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We make stringent fiducial requirements on the tau seed to ensure that the tau candidate's 
energy is well measured. The track must be fully contained in the CTC and must pass 
additional fiducial cuts similar to those imposed on the electron candidate to ensure that 
it is not incident on an uninstrumented portion of the calorimeter. Also, to suppress fake 
track contamination, we require 0.5 GeV in the tower to which the track points. The seed 
track also must be within 5 cm of the same primary vertex (|^vtx ~ z r h \ < 5.0 cm) as the 
electron track. 

In the neighborhood of the tau seed, two isolation regions are considered sepa- 
rately with different requirements made in each region: the AR < 0.175 cone and the 
0.175 < AR < 0.524 annulus. In either isolation region, a track is a shoulder track if it is 
a good quality track with > 1 GeV/c. The seed track is included in the track counting 
in the R < 0.175 cone. 

For a true hadronically decaying tau, the number of tracks in the AR < 0.175 cone 
(^conc) i s usually 1 or 3, so we require N^ e < 4. We additionally require the sum of the 
charges of the tracks in the cone to be ±1, and opposite to the charge of the electron. We 
expect the number of tracks in the annulus 0.175 < R < 0.524 around the tau seed (N^) 
to vanish in signal events, so we require jV^jf = 0. These tracking isolation cuts retain 
approximately 80% of signal events and reject 70-80% of QCD jets. 

The CES clustering algorithm used is the same as the one used in previous CDF analy- 
ses with some modifications that improve tau purity and fake rejection. In particular, 
CES clusters were formed at a larger distance in r — <p from the seed track so that this in- 
formation may be used for fake rejection. Also, a \ 2 requirement that was used in previous 
CDF analyses to ensure the consistency of the CES cluster profile with electron test beam 
data was removed here to improve efficiency without a significant sacrifice in purity. 

The algorithm forms CES clusters in the AR < 0.6 cone around the seed track by taking 
the highest energy strips (wires) in each calorimeter tower, in descending order, calling them 
seeds, and merging them with their nearest neighbors to form clusters 4 strips (6 wires) 
wide. To be a seed for a cluster, a strip (wire) must show a pulse height that surpasses 0.4 
(0.5) GeV. The cluster position is defined as the position of the center strip or wire in the 
cluster. Pulse heights from strips and wires were corrected for 77— and 4>— dependent effects, 
measured from test-beam data. The energy of a CES cluster in a tower is the CEM energy 
of that tower, weighted by the energy deposited in the CES by that cluster compared to the 
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energy deposited by all CES clusters in the tower. The predicted response in the CEM for 
a charged pion is subtracted from the energy in the tower impacted by the seed track. Wire 
clusters in the r — <p view are matched with strip clusters in the z view with similar energies 
and merged into new clusters. The cluster position must not be consistent with coming 
from the seed track. It is rejected if \rj^ ecd — ^ luster | < 0.03 and either |0 secd — clustcr | < 0.01 
or \(fi tower center — 0ciustcr| < o.Ol (the latter requirement is because the cluster position is 
assigned to the center of the tower in <p when no wire information is available). 

A CES cluster must satisfy Ex > 1 GeV to be counted as a shoulder cluster in either 
of the isolation regions. We call the number of CES clusters found in the isolation cone 
(annulus) N^ e (^ann) an d require N™ e < 3 and N^ n = 0. In addition, the CES cluster 
energies in the cone measured as 3-component vectors are combined with the measured 
momenta of the tracks to compute a tau mass, m T . We require m T < 2.0 GeV/c 2 . The cuts 
on N^ e , N^ n and m T give a combined efficiency for signal of approximately 95%. These 
three cuts additionally reject 30-50% of QCD jets. 

As described above, a tau candidate is found to be isolated through a measurement of 
track and CES cluster multiplicities in the neighborhood of a seed track. We compare the 
isolation variables among Z — > ee, Z — > and Z — > rr simulation MC samples, and in 
Z — > ee and Z — > fifi data control samples. In the Z — > ee and Z — > fifi samples, both 
lepton candidates in the event mimic the tau seed in this analysis. We find good agreement 
between data and simulation, and between electron and muon samples, in the efficiencies 
of the isolation cuts. The isolation efficiencies from the simulated Z — > tt MC sample also 
agrees with the data samples in the annulus around the tau seed (where no particles from 
tau decays are expected). 

To reduce the impact of Z/j — > ee background on the sensitivity, we require £ > 0.15 
where £ = E T ad /Sp^ j^. Here, E T ad is the hadronic energy of the tau jet cluster and Sp^ 
is the sum of the pt of all tracks with > 1 GeV/c within the cone AR < 0.175 centered 
on the jet direction. 

Since the decay products of a hadronically decaying tau are highly collimated, a tau is 
expected to leave a narrow cluster of energy in the calorimeter. We cut on the <fr — <f) and rj — rj 
moments of the jet cluster associated with the hadronic tau candidate which are defined as 
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FIG. 4: Efficiencies of the tau identification cuts applied to hadronic taus from Z — > tt events . 
follows: 



T 



E#T 

Ez e t -(vi-vo) 2 



V V 



(2) 



(3) 



The sum is over calorimeter towers in the jet cluster, and <po and t]o are the Et weighted 
center of the jet in the <fi and 7] directions. We require I^-^ < 0.1 and < 0.1. These 
two cuts together reject approximately 30-45% of QCD jets. 

Figure 0] shows the efficiency of the hadronic tau identification cuts from the simula- 
tion. We bin the efficiencies according to the visible transverse energy (Et) from the tau at 
the generator level. The total efficiency plateaus near 55% for Et ^> 35 GeV. For fiducial 
hadronic taus with Et > 10 GeV from Z decays, the average efficiency of the tau identifi- 
cation cuts is 40%. A — > tt events at = 100 GeV and tan /3 = 50 are also 40% efficient. 



D. Additional Requirements 

We make further cuts on the data sample to increase purity and further reject background. 
We require at most one recoil jet with Et > 15 GeV in the event. For further suppression 
of Z — > ee background, we reject any event with two electrons or one electron and one track 
that have a reconstructed invariant mass M ee between 70 and 110 GeV/c 2 , which removes 
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99% of Z — > ee events while retaining 90% of signal events. We require a separation of the 
tau candidates in the transverse plane, A<p > 1.5, where A0 is the azimuthal angle between 
r e and r/j. This is nearly 100% efficient for signal rejecting 20% of non-tau backgrounds as 
measured from a background-dominated data sample. 

To take advantage of the mass reconstruction technique, we divide the events into back- 
to-back and non-back-to-back samples. The full invariant mass of the di-tau system can be 
estimated only when the tau candidates are not back-to-back, as explained in more detail in 
Section IV E. The tau candidates are called back-to-back when sin A0 < 0.3, where A0 is 
the azimuthal angle between r e and Th- Note that sin A0 is the determinant of the system of 
equations which determine the di-tau mass, so when sin A0 pa 0, the solution is not unique. 

The missing transverse energy in the event, denoted IpT, is the opposite of the vector 
sum of the measured transverse energies of the event. For the non-back-to-back events, we 
use the magnitude and direction of Ifix to derive the di-tau mass. First, we define corrected 
]£t m the following way: 

g^orr = _ Et - PT j ~ ^e/' - ^ AEp . (4) 

towers muons jets 

The first term on the right side of the equation is a sum of the the transverse component of 
the energy deposited the calorimeter towers. The remaining terms improve the resolution 
by accounting for the momentum carried away by muons, energy corrections applied to the 
electron candidate, and jet energy corrections. 

It is only necessary for the simulation to correctly model T corr well in the region pp h,z > 
15 GeV/c, since that is where the mass reconstruction is utilized. We confirm that the ]fi T corr 
variable from the data is well modeled by the simulation using a sample of Z — > ee events 
with pf > 15 GeV/c. 

E. Di-tau Mass Reconstruction 

Signal events contain neutrinos that escape CDF undetected. At hadron colliders, the 
resulting energy imbalance may only be determined in the transverse plane because the 
^-component of the total momentum of the interaction is unknown. 

Nonetheless, the energy of the neutrinos from each tau decay, and thus the full mass of 
the di-tau system, may be deduced if (i) the tau candidates are not back-to-back in the 
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The contributions to the total missing energy from the leptonic and hadronic decays, 
denoted Jf,\ and Jfih , are the solution to a system of two equations and two unknowns: 



#1 sin 8i cos (pi + lp h sin 6 h cos <p h = {JP 



(5) 



# sin 0, sin 0, + $ h sin 6 h sin <f> h = ) y (6) 

Here, Jp meas is the missing energy measured for the event and 8^h are the polar angles of 
the taus and 0/^ are the azimuthal angles of the taus. The tau candidate directions are 
measured from the visible decay products. 

The reason for the first of the two criteria for the mass technique outlined above is that 
when the tau candidates are back-to-back in the transverse plane, the reconstructed mass 
is not a good separating variable because there are many high mass solutions. 

Some events may give negative solutions for Ifii and Ifih ■ We require Ifii^ > for the non- 
back-to-back events, which reduces non-di-tau backgrounds in this region. Figure El shows 
Ifii from a background-dominated sample of the data compared to simulated A°/h° — > tt 
events. The analogous distributions for Efh are similar. The Ifii^ > requirement is 55- 
60% efficient for signal, increasing with mass, while reducing non-di-tau background by 
approximately a factor of 10. When the di-tau mass is reconstructed as described next, this 
cut also improves the mass resolution. 
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FIG. 6: Di-tau mass distribution as modeled by Pythia 6.203 and the CDF detector simulation 
with parameters mA = 100 GeV and tan (3 = 50. A sin(A</>) > 0.3 cut has been imposed. 
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FIG. 7: Corrected missing energy as measured from a Z — > ee data control sample, only for events 
where > 15 GeV/c. 

We calculate the total reconstructed mass of the di-tau system using 

m2 TT = m Z/h = (Pl+Ph) 2 

= 2m 2 T + 2(# + ^r)(^ + ^D(l " cos^) 

where and ph are the 4-momenta of each tau, and E\ and Eh represent the total energy 
of each tau. E™ s and E™ s represent the energy left by their visible decay products. The tp 
is the 3-dimensional angle between the two taus. The missing energies ^ and Ifih are the 
solutions to Eqs. El and El 

Figure El shows the di-tau mass distribution reconstructed from signal events for the 
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Jet > 10 GeV/c 2 
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^> 1 Tan SppH 


1265 




> 1.5 


1117 






sin < 0.3 sin 


> 0.3 


Olli 


51 n 
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■yy- -L i Cuervo 


1 89 
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14U 


Elect Reiert 


98 


93 


m < 2 OeV/r- 2 


t7 (J 


84 


# CES 


80 


72 


Jet Width 


64 


54 


IQ| = 1 


48 


39 


Opp. Sign 


39 


28 


$i > 0, $ h > 


NA 


8 



TABLE III: A summary of the cuts imposed and the number of events remaining in the data 
sample after each successive cut. 

parameter point m^o = 100 GeV/c 2 and tan/3 = 50, for which the A or h° particle has 
an inherent width of 5.7 GeV/c 2 . The contribution from the calorimeter resolution may be 
qualitatively seen from Fig. since Z — > ee events would have a small Jfir- The remaining 
contribution to the width of the di-tau mass distribution comes from the approximations 
that are needed to implement the technique, listed above. 

Table III summarizes the cuts that we impose and the number of events remaining in the 
sample after each cut. Table IV summarizes the efficiency of each cut on the Z — > tt and 
A + h° — > tt simulated events. 
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Cut 



Efficiency (%) 





Z — > rr 


A° + h° -» rr 


BR(r -> e) 


32.5 




OZ.O 


|»7e| < 1-2 


19.8 




ZD.O 


Electron ID & AT iso 


= 35.6 




OI.O 


L2 Trigger 


90.8 




qi n 
yi.u 


Z Rejection 


93.3 




on q 

yu.o 


N jets < 3 


97.7 




yo.4 


Fiducial Jet 


35.3 






Jet E T > 10 GeV 


89.4 




on % 

yu.o 


> 1 Tau Seed 


64.0 




UtJ. ± 


A4> > 1.5 


100.0 




1 nn n 




sin Acf) < 0.3 


sin 4 


\/h ^ n ^ 

_A(J9 ^> U.O 


sin Ac/) 


84.1 15.9 


sn 7 

oU. ( 


1 Q /I 

ry .4 


# Tracks 


84.3 85.9 


oZ.o 


OO.U 


Elect. Reject 


95.5 96.0 


t)o.Z 


yo.u 


m r < 2 GeV 


98.5 98.6 


98.5 


98.4 


# CES 


97.0 96.6 


96.9 


97.7 


Jet Width 


92.2 92.5 


92.7 


92.0 


|QI = 1 


93.1 92.9 


92.7 


92.8 


Opp. Sign 


99.8 100.0 


99.8 


100.0 


$i > 0, $ h > 


NA 56.8 


NA 


58.0 



TABLE IV: Efficiency of each cut imposed on Z — > rr and A + h° — ► rr simulated events. Back- 
to-back and non-back-to-back events are denoted sinA0 < 0.3 and sinA</> > 0.3, respectively. 



V. BACKGROUNDS 

The backgrounds can be classified into two categories: physics and instrumental back- 
grounds. The former includes Z — > rr and Z/'-f* — > ee. The latter fall into three categories: 
events containing a conversion pair and a recoil jet, W(— > ez/)+jets events and events con- 
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FIG. 8: Hadronic tau fake rates as measured in the lepton control samples compared with that 
measured in the jet control samples. 

taining a jet that fakes an electron. All of these contain a fake hadronic tau. We model signal 



and physics backgrounds using Pythia Monte-Carlo and the CDF detector simulation |21|. 

We estimate the rate of all fake tau contributions combined instead of estimating each 
source of fakes separately. This way, we do not rely on modeling of jets, nor on limited 
control samples for each separate background. We estimate the number of fakes expected 
to pass all of the cuts using the fake rate technique described next. 

We use five different control samples from CDF data to measure the rate at which a 
jet fakes a hadronic tau by passing the tau identification cuts for this analysis. They are: 
(i) a sample dominated by events where a photon produced an electron pair, (ii) a sample 
dominated by W(— > e^)+jets, (iii) another dominated by W(— > yUz/)+jets and two samples 
of events collected using inclusive triggers with a jet satisfying (iv) > 20 GeV and (v) 
Et > 50 GeV. The first three are classified as lepton samples and the latter two are jet 
samples. 

A fake rate is defined as the probability that a jet passes the hadronic tau identification 
requirements as described in Sec. IIVCI For all data control samples, we measure the fake 
rates in bins defined by the Et of the jet associated with the tau candidate. Figure |Hl shows 
the fake rates measured from all three lepton samples combined and from both jet samples 
combined. We measure fake rates less than about 1% from the jet samples, which are as 
good as in previous tau analyses J3] • We find that the fake rates measured from the lepton 
samples are approximately a factor of 2 higher than from the jet samples. Two of the lepton 
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samples are dominated by W+jets events where the recoil jet comes from a quark; quark jets 
are narrower, have lower multiplicity and are thus more likely to pass the tau identification 
cuts [37]. Still, the reason for this difference is not definitively understood. Since the analysis 
is performed in a lepton sample, we use the fake rates measured from the lepton samples 
as the central value and the difference between the rates from the two different types of 
samples as a systematic uncertainty. The histogram in Fig. |H] shows the Et distribution of 
jets in the data sample just before tau identification cuts are applied. The fake rates are 
folded into this Et distribution to predict the rate of fake tau background. 

We find that fake taus from the 1^+jets samples are opposite-sign from the lepton in 
the event (67.1 ± 3.0)% of the time. This is due to a correlation between the W and the 
recoiling quark. The isolation cut enhances this correlation due to charge conservation. In 
the conversion sample, the opposite-sign requirement is consistent with being 50% efficient. 
We take the central value of the efficiency for the opposite-sign cut to be the average of the 
two ((67.1 + 50)/2 = 58.6%) and equally likely to be anywhere between 50.0% and 67.1%. 

To estimate the number of fake taus expected to pass all of the cuts, we use the data 
sample with the following cuts removed: (i) tau ID cuts (ii) the opposite sign requirement 
and (iii) $y 1)T2 > 0. Then, we fold in the measured fake rates with the Et spectrum of jets. 
There are 6478 events that pass these cuts and which contain at least one jet passing the 
fiducial cuts. We apply our measured fake rates to 6972 jets from these events. For each 
jet, we give it a weight equal to the fake rate measured for its Et- We expect 21.0 ± 12.1 
back-to-back and 29.8 ± 16.8 non-back-to-back for a total of 50.8 ± 29.0 fakes before the 
remaining cuts are applied. 

We apply two final cuts to improve the purity of the sample. We apply the opposite 
sign requirement, taking the efficiency to be 58.6%. Including the systematic error from 
this cut, we expect 10.5 ± 6.0 back-to-back and 14.9 ± 8.5 non-back-to-back events at this 
stage. The final cut is the Jfir ltT2 > cut applied to the non-back-to-back events only. This 
is measured from the same background-dominated sample, subtracting out Z — > rr and 
Z — > ee contamination. We find that this cut removes 89.0% of the fake tau background. 
This brings the number of predicted non-back-to-back events with a fake tau to 1.6. Added 
to the back-to-back events, we expect 1.6 + 10.5 = 12.1 events containing a fake hadronic 
tau to pass the analysis cuts. At each stage of this estimation, we account for the 10%-level 
Z — > rr and Z — > ee contamination of the background-dominated sample. 
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VI. SUMMARY OF SYSTEMATIC UNCERTAINTIES 



In Tab. V , we summarize the systematics on the backgrounds and signal. Z — > ee is not 
included in the table because the expected rate is based on a low number of background 
events. We expect 0.6 ± 0.3 Z — > ee events in the counting experiment. The systematics 
on the fake tau background are described in Sec. 

3q . The systematic error on the 



The error on the Run lb luminosity at CDF is 4.1% 
electron identification cuts is 1.8%, as noted in Sec. IIV1 This error comes from the limited 
size of the Z —>■ ee sample used for comparing the simulation with data. The systematic 
error due to the modeling of the trigger efficiency is obtained by moving the parameters 
in the energy and //-dependent efficiency by one standard deviation in each direction and 
measuring the effect on the total efficiency of all cuts. This systematic is 1.6% (1.7%) for 
A°/h° ^tt (Z^tt). 

The one systematic uncertainty on the yield of signal events that varies significantly with 
the mass of the Higgs boson is the uncertainty on the cross section due to the low a/Q 2 
tail. Here are the systematic uncertainties due to this effect for each parameter space point 
considered in this note: m^o = 100 GeV/c 2 , tan/3 = 50, 0.5%; m^o = 110 GeV/c 2 , tan/3 = 
50, 2.5%; m A o = 120 GeV/c 2 , tan/5 = 50, 3.5%, m A o = 140 GeV/c 2 , tan/5 = 50, 7.2%; 
m A o = 140 GeV/c 2 , tan/5 = 80, 21.3%. The systematic uncertainty on the Z production 
cross section is 1.7% based on the CDF measurement 3^]. 

By studying isolated pions, we have shown that the detector simulation overestimates 
the width of jets from charged pions, and thus underestimates the efficiency of the cut on 
jet width on the hadronic tau candidates. Therefore, we take the efficiency for that cut to 
be the average of the simulated efficiency and 100%, which is (100. 0+92. 2)/2 = 96.1%. The 
systematic error on this cut is half the difference between the efficiency from simulation and 
100%, which is (100. 0-92. 2)/2 = 3.9%. The electron rejection cut applied to the hadronic 
tau candidate is not modeled sufficiently for hadronic decays. The same study of isolated 
pions previously mentioned showed that the hadronic energy deposited by charged pions 
is underestimated by the simulation, so that the efficiency of this cut is underestimated. 
We take the efficiency to be the average of the efficiency from simulation for hadronic taus 
and 100%, which is (100. 0+95. 5)/2 = 97.8%. The systematic error on this cut is half the 
difference between the simulated efficiency and 100%, which is (100. 0-95. 5)/2 = 2.2%. We 
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systematic uncertainty 


A°/h° 


■ rr Z — > rr 


Fakes 




(%) 


(%) 


(%) 


luminosity 


4.1 


4.1 




cross section 


0.5 


1.7 




electron ID 


1.8 


1.8 




sample dependence of fake rates 






57.1 


opposite sign 






14.7 


jet width 


3.9 


3.9 




jet energy scale 


1.0 


1.2 


_ 


trigger efficiency 


1.6 


1.7 




electron rejection 


2.2 


2.2 




total error (%) 


6.7 


7.0 


NA 



TABLE V: Summary of systematic uncertainties on the counting experiment for the hnal cuts. 
Here we refer to Higgs boson events corresponding to the parameter space point m^o = 100 GeV, 
tan/3 = 50. All systematic uncertainties are quoted as a percentage of the total number of events 
observed after all cuts are applied. Z — > ee is not included in the table because it is based on a 
low number of MC events. The errors attributed to the rate of fake taus are not both Gaussian, 
and are therefore not added in quadrature. See the text. 

obtain the jet energy scale systematic by varying the energies of all of the jets (except those 
identified as electrons) up and down by 5%. The resulting systematic error is 1.0% (1.2%) 
for A°/h° — ► tt (Z — > rr). There is a systematic uncertainty attributed to the tau fake 
rates since we measured different rates in lepton and jet samples. We set this systematic 
uncertainty to the the difference between the fake rates measured in the two types of samples; 
it is the dominant systematic error at 57.1%. 

VII. RESULTS 

We set limits on direct A°/h° production in the MSSM based on a counting experiment 
using events from both the back-to-back and non-back-to-back samples. Then we show the 
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Number of Tracks 



FIG. 9: Track multiplicity distribution for the counting experiment, but without the \Q\ = 1, 
opposite sign and N^fe < 4 requirements. 

limits achieved from a binned likelihood fit to the di-tau mass distribution from the non- 
back-to-back events alone. Our nominal limits come from the counting experiment utilizing 
both back-to-back and non-back-to-back events. 

A. All Events 

We plot the track multiplicity of the tau candidates after imposing all of the analysis cuts 
except the following cuts: i) | ^ Q%\ — 1, ii) < 4 and iii) the opposite sign requirement. 

Figure El shows the number of tracks in the 0.175 cone around the tau seed in the hadronic 
tau candidate in the event. We expect 78.2 events to appear in this plot and we observe 
81. Table IVT1 lists the number of events expected of each background type. The data and 
the prediction show good agreement at each stage. Of the 47 final observed events, 35 are 
1-track and 12 are 3-track. The final three cuts reduce the fake background in these two bins 
by nearly a factor of 2 compared to that shown in Fig. Eland leaves the Z — > tt background 
in those bins virtually unchanged. 

B. Non-Back-to-Back Events 

Figure El shows the track multiplicity distribution for the non-back-to-back events only 
(since we will be performing the di-tau mass reconstruction on these events). These events 
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Z — ► tt Z — > ee Fake Taus Total Obs. 



All cuts but 

\Q\ = 1, Opposite Sign 46.1 0.56 31.5 78.2 81 

\Q\ = i,# Tracks< 4 42.4 0.56 20.2 63.1 61 

Opposite Sign 42.4 0.56 11.8 54.8 47 

TABLE VI: Number of each background type expected, in the track multiplicity plot and after 
subsequent cuts are imposed, compared with the number observed. 
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FIG. 10: Track multiplicity distribution for the non-back-to-back events with the Jfiij%> require- 
ment, but without the < 4, \Q\ = 1 and opposite sign requirements. 

have passed the $i,h> requirement but the \Q\ = 1, opposite sign, and < 4 require- 

ments have not been imposed. We expect 9.2 events to appear in the plot and we observe 
15. The probability for the disagreement between data and the prediction to be equal to 
or more than what was observed is 3.2%. Table IVTT1 lists the number of events expected of 
each background type in Fig. El including non-back-to-back events only. We also show the 
expectation compared to the observed after the subsequent cuts are imposed. The final 8 
events contain 5 1-track events and 3 3-track events. 
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Z — > TT 


Z — > 66 


Fake Taus Total Obs. 




All cuts but 








\Q\- 


= 1, Opposite Sign 4.4 


0.05 


4.5 


9.2 15 


\Q\ 


= 1,# Tracks< 4 4.1 


0.05 


3.0 


7.1 13 




Opposite Sign 4.1 


0.05 


1.7 


5.9 8 



TABLE VII: Number of each background type expected, in the track multiplicity plot and after 
subsequent cuts are imposed, compared with the number observed. Only non-back-to-back events 
are included. 

C. Limits 

The numbers of observed events do not show an excess above the standard model expec- 
tation for background processes. Therefore, we set a limit on the product of the cross section 
and branching ratio of A°/h° — > tt (cr- BR) at 95% confidence level. We take the branching 
ratio to be 9% for all parameter space points considered. We use a Bayesian method with 
a flat prior based on a likelihood that is smeared to account for systematic errors. 

Table IVTlTl shows the predicted and observed upper limits on a- BR for A°/h° production 
in the MSSM as a function of m^o at tan/5 = 50. For m^o = 100 GeV/c 2 , tan/5 = 50 in the 
MSSM, the expected limit is 23.4 signal events, corresponding to 102 pb. Since we observed 
slightly fewer events than we expected, the observed limits are better than our expected 
limits. The observed limit is 77.9 pb. The limits on the a- BR improve with increasing mass 
since the efficiency improves, but we are less sensitive to the MSSM theory at higher mass 
due to the steeply falling predicted cross section. 

The cross section for producing A°/h° in the MSSM scales with (tan/3) 2 . The sen- 
sitivity does not improve by that same factor, however, because the Higgs boson width 
scales as (tan/5) 2 , and the tail at low also becomes more prominent with increasing 
tan /5, increasing the systematic error due to the uncertainty in the cross section in this 
region. At m A o = 140 GeV/c 2 and tan/3 = 80, the uncertainty on the efficiency of the 
selection of Higgs boson events due to this low mass tail is 20%, compared to 7.2% at 
77^40 = 140 GeV/c 2 , tan/5 = 50. Also, both effects bring down the efficiency at higher tan/5: 
at m^o = 140 GeV/c 2 and tan/5 = 80, the efficiency is similar to the efficiency at a lower 
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Mass 


cr-BR 


Eff. 


95% CL 


95% CL 


GeV/c 2 


pb 


% 


Exp. (pb) Obs. (pb) 


100 


11.0 


0.78 


102 


77.9 


110 


5.8 


0.91 


87.4 


67.2 


120 


4.0 


0.97 


82.3 


63.0 


140 


1.6 


1.1 


75.2 


57.9 



TABLE VIII: 95% CL limits vs. m^o for tan (3 = 50. The first two columns are the mass of the 
Higgs boson for each parameter point and the corresponding theoretically predicted a- BR. 



tan (3 


a- BR Eff. 


95% CL 95%CL 




pb % 


Exp. (pb) Obs. (pb) 


50 


17.8 1.1 


75.2 57.9 


80 


65.9 0.85 


118 90 



TABLE IX: 95% CL limits vs. tan (3 for m A o = 140 GeV/c 2 . All limits are quoted in pb. The 
efficiencies shown do not include branching ratios. The first two columns are the mass of the Higgs 
boson for each parameter point and the corresponding theoretically predicted a- BR. 

mass point: m^o = 110 GeV/c 2 , tan/3 = 50. Table HX1 shows the limits for two different 
values of tan/3 for m^o = 140 GeV/c 2 . These limits are also summarized in Figure ITTI 

In the non-back-to-back region, after all cuts are applied, we expect 5.9 events and observe 
8. Figure ^] shows the di-tau mass distribution for these 8 events compared with the 
expectation. In this region, limits are obtained by fitting the mass distribution using a 
binned likelihood, with 14 bins between 60 GeV/c 2 and 200 GeV/c 2 in di-tau mass. 

Table X shows the upper limits on the MSSM o- BR obtained from the binned likelihood 
for four different values of m^o for tan (3 = 50. At m^o = 100 GeV/ c 2 , where the Higgs boson 
mass is nearly on top of the Z mass, the mass reconstruction is less effective than at higher 
masses, so the expected limit from the binned likelihood in the non-back-to-back region is 
approximately 2.4 times worse than the expected limit from the counting experiment without 
the non-back-to-back requirement. At m^o = 140 GeV/c 2 , which is approximately 2 RMS 
away from the Z in the di-tau mass variable, the expected limit from the binned likelihood 
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FIG. 11: Summary of limits achieved in this analysis on the cross section for directly produced 
Higgs bosons for each parameter space point considered. 
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FIG. 12: Di-tau mass reconstruction for events in the non-back-to-back region. It is compared 
with the background prediction, shown for fake tau backgrounds and then the Z — > tt background 
added to it. The negligible contribution from Z — > ee is not shown. 

using the non-back-to-back events is 2.1 times worse than the counting experiment limit, 
showing a modest improvement, but still not coming close to the expected limit from the 
counting experiment. With more data collected in Run 2, the power of the di-tau mass 
reconstruction technique will improve. Table XI shows the upper limits obtained from the 
binned likelihood for two different values of tan (3. 
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Mass 


a- BR 


Eff. 


95% CL 


95% CL 


GeV/c 2 


pb 


% 


Exp. (pb) Obs. (pb) 


100 


11.0 


0.093 


247 


395 


110 


5.8 


0.10 


218 


379 


120 


4.0 


0.105 


199 


363 


140 


1.6 


0.12 


158 


309 



TABLE X: 95% CL limits vs. m^o for tan j3 = 50, using the binned likelihood. The first 
two columns are the mass of the Higgs boson for each parameter point and the corresponding 
theoretically predicted a- BR. 



tan/? a- BR Eff. 95% CL 95% CL 
pb (%) Exp. (pb) Obs. (pb) 
50 1.6 0.12 158 309 
80 5.9 0.09 274 533 

TABLE XI: 95% CL limits vs. tan (3 for m A o = 140 GeV/c 2 , using the binned likelihood. All limits 
are quoted in pb. The efficiencies shown do not include branching ratios. The first two columns are 
the mass of the Higgs boson for each parameter point and the corresponding theoretically predicted 
a- BR. 

VIII. CONCLUSIONS 

We have performed a search for directly produced Higgs bosons decaying to two taus 
where one tau decays to an electron and the other hadronically in Run lb data at CDF. 
This is the first Higgs boson search based on the di-tau final state at a hadron collider. The 
number of events that pass all of the cuts is consistent with the background expectation. This 
agreement between data and background demonstrates our capability to reconstruct Z — > rr 
final states, the irreducible background for this analysis. At a benchmark parameter space 
point, rrijsfi = 100 GeV/c 2 and tan/3 = 50, we are sensitive to a a- BR of 102 pb compared 
to the 11.0 pb predicted in the MSSM. The observed limit at this parameter space point is 
77.9 pb. 
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A di-tau mass reconstruction is performed for tau candidate pairs which are not back- 
to-back, for the first time with hadron collider data. The modest sensitivity that one gains 
from a limit binned in mass is not nearly enough to make up for the hit in efficiency taken 
when only non-back-to-back events are considered. At m^o = 100 GeV/c 2 , the binned mass 
limit at tan/3 = 50 from non-back-to-back events alone is 395 pb. At m^o = 140 GeV/c 2 , 
tan p = 50, the binned mass limit is 309 pb, showing a modest improvement as the Higgs 
boson mass is moved away from the Z mass. 

While this search does not have the sensitivity to the Standard Model Higgs boson that 
prior searches using decays to pairs of b quarks, it lays the groundwork for similar analysis 
to be performed by future experiments. The di-tau mass reconstruction technique demon- 
strated here may also be useful for searches for other processes where a Higgs boson is 
produced with a recoil, such as Hb or Hbb. 
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